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Abstract 

We discuss the finite volume corrections to fx and Bk by using the one-loop chiral 
perturbation theory in full, quenched, and partially quenched QCD. We show that the 
finite volume corrections to these quantities dominate the physical (infinite volume) 
chiral logarithms. 



PACS: 12.39. Fe (Chiral Lagrangians), 11. 15. Ha (Lattice gauge theory) 



1 Introduction 



Because of the limited computing power, current lattice computations of the hadronic ma- 
trix elements involving kaons are plagued by the necessity for introducing three important 
approximations: 

1. The (partially) quenched approximation. 

2. The extrapolation in the light quark masses: because of the inability to simulate di- 
rectly with the physical u/d quarks, one works with masses not lighter than about one 
half of that of the physical strange quark and then extrapolates to the physical niu/d- 
Once the lighter quark masses are probed on the lattice, such as those announced in 
ref. [1], the finiteness of the lattice volume becomes also an important approximation. 

3. Degeneracy of valence quark masses in the kaon: matrix elements involving kaons are 
obtained with "kaons" consisting of degenerate valence quarks whose mass is tuned 
in such a way as to produce a pseudoscalar meson with its mass equal to the physical 
rriKO = 0.498 GeV. 

In view of the great importance of the — mixing amplitude in constraining the shape 
of the CKM unitarity triangle [2], a quantitative estimate of the systematic errors induced 
by the above listed approximations is mandatory. That is where chiral perturbation theory 
(ChPT) enters the picture and offers a systematic approach for quantifying (at least roughly) 
the size of these errors. In ChPT one computes the coefficients of the chiral logarithms for 
various hadronic quantities in order to (a) examine whether or not the quenched approxima- 
tion introduces potentially large systematic error, (b) guide the chiral extrapolations, and (c) 
quantify the impact of the degeneracy in the quark masses on the evaluation of the hadronic 
matrix elements. The coefficients of the chiral logarithms are predicted by quenched and full 
ChPT. Although convincing evidence for the presence of chiral logarithms in any numerical 
lattice data is still missing, a slight discrepancy from the linear (or quadratic) dependence 
on the variation of the light quark mass is occasionally observed. Before identifying such 
a discrepancy as an indication of the presence of chiral logarithmic behavior, one should 
make sure that the effects of finite volume are well under control. In particular, we would 
like to know how the finiteness of the lattice volume modifies the chiral logarithmic behav- 
ior of fx and B^- In this paper we present the expressions obtained in three versions of 
ChPT that are relevant to present and future lattice simulations, i.e., in quenched ChPT 
(QChPT), partially quenched ChPT (PQChPT) and in full (standard) ChPT. Concerning 
the PQChPT, we will consider the case of A^sea = 2 degenerate dynamical quarks, which 
is the current practice in the lattice community. Those expressions, obtained in both the 
finite and infinite volumes, are then used to (i) show that chiral logarithmic behavior of fx 
and Bk is indeed modified by the finiteness of the volume, and (ii) assess the amount of 
systematic uncertainty induced by the finiteness of the lattice volume. As expected, finite 
volume effects increase as the mass of the valence light quark in the kaon, m^, becomes 
smaller (we keep the strange quark mass fixed to its physical value). For quark masses 
iriq > ms/3 and volumes V ^ {2 fm)'^, the finite volume effects are negligible. We will argue 
that even if one manages to push the quark masses closer to rriu/di finite volume effects will 
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start overwhelming the effects of the physical (infinite volume) chiral logs (unless one uses 
very large volumes). This unfortunately complicates the efforts currently being made in the 
lattice community to observe the chiral log behavior directly from the lattice data. Our 
finite volume ChPT formulae for fx and Bk may be used to disentangle the finite volume 
effect from the physical chiral logarithmic dependence. Obviously this can be done only if 
the volume is sufficiently large and thus the finite volume corrections safely small enough 
to justify neglect of the unknown higher order corrections in the chiral expansion. 

The remainder of this paper is organised as follows. In Sec. 2 we compute the chiral log 
corrections to fx and Bk in all three versions of ChPT in infinite volume. In Sec. 3 we 
discuss the same chiral corrections but in the finite volume. Both these sets of expressions 
are then combined in Sec. 4 to examine the impact of the finite volume artefacts on the 
chiral behavior of fx and Bx- In Sec. 5 we discuss the finite volume effects on fx and Bx 
and we briefly summarise in Sec. 6. 



2 Results in infinite volume 

To simplify the presentation and for an easier comparison with results available in the 
literature, we first briefly explain the notation adopted in this work and then present the 
expressions for the chiral logarithmic corrections that we computed in all three versions of 
ChPT. 



2.1 Chiral Lagrangians 

For the full (unquenched) ChPT we use the standard Lagrangian [3,4]: 

/:chPT = y tr [(9^St) (9^S) + + X^S ] , 
with / being the chiral limit of the pion decay constant, /tt = 132 MeV. In addition. 
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For the calculation in QChPT we will use the Lagrangian introduced in refs. [5, 6]: 
>CQChPT = ^str [(9,,St)(9^S) + Stx + x^S] -ml^l + ao{d^%){d>'%) . 



(1) 



(2) 
(3) 



(4) 
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where $o = str[$]/A/6, is proportional to the graded extension of the r]', the trace over the 
chiral group indices has been replaced by the supertrace over the indices of the graded group 
SU{3\3)i X SU{3\3)ji, and the fields S and x ^ow graded extensions of S and Xy just 
defined above. 

Finally, we choose the SU{5\3)l x SU{5\3)r setup for the PQChPT, i.e., three valence 



quarks (m, d, s), with masses rriq 



m,, 



rud 7^ ms, and two degenerate sea quarks (wsea, c?sea) 



of mass msca- The Lagrangian is of the same form as the quenched one in eq. (4), except 
that the indices now run over the graded group SU{5\3)l x SU{5\3)r, and the fields S and 
X are extended to include the sea-quark sector [7]. Moreover, because of the presence of sea 
quarks, the rj' decouples and $o can be integrated out of the Lagrangian [8]. 

Throughout the paper, the evaluation of the chiral loop integrals is made by using naive 
dimensional regularisation and the so called "MS + 1" renormalisation scheme of ref. [3]. 



2.2 One-loop chiral log corrections to /a and B 



K 



We begin by collecting the ChPT expressions for fx and Bk in infinite volume. We adopt 
the standard definition of the Bk parameter, namely. 



B 



K 



(ir°|g7/.(l-75)rf^7M(l-75)c?|i^°) 



(5) 



which is equal to 1 in the vacuum saturation approximation. The bosonised version of the 
relevant left-left (AS" = 2) operator reads 



927 (S9,St)^^ (S9''St)^ 



(6) 



To compute the chiral loop corrections to fx, we use the standard bosonised left handed 
current: 



J. 



s7m(1 -75)c? 



4 (^^'^^') 



ds 



(7) 



In the following we will leave out the analytic terms (those accompanied by low energy 
constants) and focus only on the nonanalytic ones. As we will see, the analytic terms are 
not relevant to the discussion of finite volume effects. 
The chiral logarithmic corrections to fx are 
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where ChPT, PQChPT, and QChPT stand for the full, partially quenched (iVg, 
quenched chiral perturbation theory. In the above formulae, 



m 



ss 



2Boms 



m, 



22 



2BQms = 2m 



K 



m^ 



m 



23 ^ 5o(m^ + msca), m^g = 5o(mg + msea) 



(10) 
2), and 

(11) 



We stress that we work in the exact isospin symmetry limit, i.e., = = m^- The 
results listed above agree with the ones available in the literature: eq. (8) was first obtained 
in ref. [3], eq. (9) in refs. [9,10], and eq. (10) in refs. [5,6]. 
For the Bk parameter, we obtain: 



B 



K 



rstrec 



ChPT 



(47r/)2 
[Tmj^ — ml 



+ log 

2 



m 



K 



mt (ml- + ml) , / m 



4 m|- 



log 



m: 



4m|- 



-log 



(12) 



B 



K 



N PQChPT 



rjtree i 



1 - 



m 



ss 



m„ 



m\ + ml 
2m\ 



m 



-— m 
2 



K 



SS' 



2m\ 



^ + ml^ ^ 



m 



K 



mt 



m 



K 



log 



log 



m 



22 



mt 



m 



K 



2 log 



22 



(13) 



B 



K 



QChPT 



Rtree 



3(47r/)^ 



6m ^ + 6m ^ log 



m 



K 



+3 m: 



m22 + m 



22 



m 



log 



m 



22 



mn 



+ 3 m: 



m\ + ml_ / m|' 



m 



K 



m|^ + m^r, m: 



22 



mj^(m 



K 



mt 



log 



m: 



22 



mf 



-2ao 



3m^ — 



ml 



+ 



m; 



m]^ + m 



2 4 

ml — m„ 



m 



K 



m 



K 



m 



K 



mt 



-log 



ml 



m 



22 



m]^ + mT^- m,, — m: 



22 



'■22 



m 



"^1- - "^22 



log 



m. 



22 



-4 



(14) 



These results agree also with the ones previously computed in full ChPT [11, 12], PQChPT [10], 
and QChPT [6, 13], where more details about the actual calculation can be found. 



3 Results in finite volume 

The calculation of the chiral logarithmic corrections in a finite box of volume V = Lr", with 
periodic boundary conditions, is completely analogous to that in infinite volume, except for 
the fact that loop integrals now become sums over discretised three-momenta. As on the 
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lattice, at the end of the calculation, the times of the kaon fields in the correlation function 
are sent to infinity. To abbreviate the expressions, we first introduce 

2 ^2,2 2 ^2,2 2 ^2,2 n tx\ 

uj^ = q + , ujj^ = q + ruj^ , ^22 = 5 +"^22^ (15) 

and analogously for tuis, UJ23, with the corresponding masses already defined in eq. (11). As 
in infinite volume, mo and ao are the t]' parameters of the quenched theory. 
For the decay constant J'k in all three versions of the ChPT, we obtain. 
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while for the Bx parameter we have 
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now faced with the problem of evaluating the sums over discrete momenta 
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3.1 Evaluation of the chiral loop sums 

The sums that appear in the calculation of the tadpole diagrams are of the form 
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where M stands for the generic mass. It is very easy to verify that 
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For finite L, one can write 
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where C,s{L,M) is simply the difference between the finite volume sum and the infinite 
volume integral. This function is finite and needs no regularization since it represents 
an infrared effect. In other words, the integral and the sum diverge in the same way. 
Equation 24 can then be considered as a way to regularize the sums which, in addition, 
allows us to adopt the same renormalization scheme for both integrals and sums. 
In the following few steps, we show how ^s(L, M) is simplified to 
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where the eUiptic theta function 'i9(r) is defined as ^ 

oo 
n=— oo 

and satisfies the Poisson summation formula [15] 



Applying the formula (27) to eq. (25), we get 
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In the asymptotic limit L oo, the theta function behaves as i?(L^/4r) ~ 1 + 2e ^^/^^, so 
that in the same limit we can write 
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4 Impact of finite volume eff^ects on the chiral behavior 

of fx and Bx 

In recent years considerable effort has been invested in controlling the chiral extrapolations 
of the hadronic matrix elements computed on a lattice. To guide the extrapolation from the 
directly accessible quark masses, r ^ 0.5, down to the physical r — >• r^/d = 0.04, one can 
rely on the expressions obtained in ChPT (quenched, partially quenched, or full). Those 
expressions, however, contain chiral logarithmic terms which so far have not been observed 
in the numerical studies. An important task before the lattice community is to lower the 
quark mass and get closer to the region in which the chiral logarithms become clearly visible. 
However, by decreasing the quark mass, the sensitivity to the finiteness of the lattice box of 
side L becomes more pronounced. Moreover, the finite volume effects modify the nonlinear 
light quark dependence in the same way, i.e., they enhance the chiral logs. The problem is 
that the nonlinearity induced by the finite volume is larger than that due to the presence 
of physical chiral logarithms. To illustrate that statement, in fig. 1 we plot the chiral log 
contributions in the finite and infinite volumes, by using the expressions presented in the 
previous section for both and Bk in all three versions of ChPT. From that plot we see 
that it is very difficult to distinguish between physical chiral logarithms (thick curves) and 
the finite volume effect, even if one manages to work with very light quarks on the currently 

oo 

^The function '!9(r) is obtained from the eommonly used function '!?3(m, q) = g" ^2nui ^ after replac- 

n— — oo 

ing, u — and q ~ e^^ . For the numerical analysis, we use the function predefined in Mathematica, 
namely, EllipticTheta[3, 0, e^'^] . For more details on the elliptic functions, see ref. [14]. 
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Figure 1: From top to bottom, we plot the chiral logarithmic corrections as predicted in full, partially 
quenched (vsea = mgca/'m^^^^ = 0.5j and quenched ChPT, respectively, as functions of the light valence 
quark mass r = mq/ms, where the strange quark mass is fixed to its physical value. In each plot the thick 
line corresponds to the physical (infinite volume) chiral logarithm, whereas the other four curves correspond 
to the logarthmic contributions computed in the finite volume V = L^, where for L we choose the values 
shown in the legend. The renormalisation scale is chosen to be fj, — 1 GcV. 
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used lattice volumes. For smaller masses, at which the chiral logarithms are expected to set 
in, the finite volume effects completely overwhelm the physical non-linearity. 

A possible way out would be to fit the lattice data to the finite volume forms (see Sec. 3) 
and not to those of the infinite volume, given in Sec. 2. That, of course, is legitimate if 
one assumes the validity of the NLO ChPT formulae. Finally, the curves corresponding to 
L = 1 fm should be taken cautiously because this volume may be too small for ChPT to 
set in, as recently discussed in ref. [16]. 



5 Finite volume corrections 

In this section we combine the formulae derived in Sees. 2 and 3 to discuss the shift of fx 
and Bk induced by the finite volume effects. Before embarking on this issue, we first briefiy 
remind the reader about the similar shift in the case of the pion mass where, for large L, 
the one-loop ChPT expression indeed agrees with the general formula derived by Liischer 
in ref. [17]. Since the analogous general formulae for fx and Bk do not exist, we will derive 
them by taking the large L limit of our one-loop ChPT formulae. 

5.1 Contact with Liischer 's formula 

To make contact with Liischer's formula, we subtract the one-loop chiral correction to the 
pion mass (squared) as obtained in the full (unquenched) ChPT in infinite volume from 
the one obtained in a finite volume. By evaluating the sums, as described in the previous 
section, we have 
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which coincides with the result of ref. [18]. ^ Analytic terms in ml{oo) and ml{L) were 
omitted since they cancel in [Ami) /ml- 

To recover Liischer's formula, one takes the limit L ^ oo, which amounts to using the 
asymptotic form (29) in eq. (30), 
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^Notice that the function gr{M'^ ,0, L), defined in ref. [18], is related to (,s{L,M) through 
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where only the leading exponential has been kept. The benefit of eq. (30) is that it offers 
insight into the subleading terms, suppressed by higher powers in e""*'"^ in Liischer's for- 
mula. In the range of volumes in which the right-hand side of eq. (31) becomes sizable, the 
subleading exponential terms cannot be neglected and the formula (30) has to be used. For 
the volumes currently used in lattice simulations, these corrections are important if one is 
to work with very light pions. 

Before closing this subsection, two important comments are in order, though. First, 
Liischer's formula relates the finite volume shift of the pion mass to the vr — tt scattering 
amplitude. Equation 31 refers to the tree level n — ir scattering amplitude. A recent study 
in ref. [16] shows that the inclusion of the NLO chiral corrections to the tt — vr scattering 
amplitude produces a sizable correction to eq. (31). It is, however, not clear whether or 
not such a conclusion persists in the full (nonasymptotic) case, i.e., eq. (30). It is even 
less clear if such a conclusion carries over to other quantities. To resolve that issue, one 
should compute the finite volume two-loop chiral correction to the pion mass (and to other 
quantities), which is beyond the scope of the present work. It is clear, however, that before 
this point is clarified, one cannot safely use the one-loop calculation to correct for the 
finite volume effects. At present the one-loop ChPT finite volume expressions are useful 
for making a rough estimate of the finite volume corrections. The second comment is that 
the derivation of Liischer's formula relies crucially on unitarity. Since the unitarity in the 
partially quenched and quenched theories is lost, Liischer's formula is not valid in these 
theories. 



5.2 Finite volume corrections to / 
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We now use the expressions for the decay constant fx, derived in the infinite [eqs. (8)-(10)] 
and finite [eqs. (16)-(18)] volume cases, to estimate the shift in fx due to the finiteness of 
the volume. For that purpose we define 
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It should be clear that the analytic terms multiplied by the low energy constants (omitted 
in Sees. 2 and 3) cancel in the ratio (32). ^ Finally, we have 
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■^Written schematically. fK{oo) — f^^'°'^{l + log^+Cmq) is equivalent to = /^^ (oo)(l — log^^ — Cm^), 
where the analytic term is multiplied by the generic low energy constant C. is the same in finite and 



infinite volumes, so that one simply obtains fxiL) 
effects of low energy constants cancel. 



fKioo){l — logg^ + logj^). Thus, at this order, the 
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Fi gure 2: The Gnite volume corrections to fx in full, partially quenched and quenched theory, 
eqs. (33), (34), (35), respectively. The partially quenched case that we consider is the one with iVf = 2 
dynamical quarks degenerate in mass for which we take Vgca ~ msea/'mP^'y^ = 1,0.5,0.2. Each plot corre- 
sponds to a different value of the size of the side of the box L, indicated in the plots. We keep the same 
scale, to better appreciate the reduction of the finite volume effects as L is increased. 



■.r N QChPT . fn, rr,2 



,2 



2 ^1 ^22) - " / " ^1 m.) I . (35) 

Similar expressions in full and quenched ChPT, but for the case in which the kaon consists 
of two quarks degenerate in mass, were obtained in ref. [6]. For the general nondegenerate 
case and for PQChPT, the above formulae are new. In ref. [19], the finite volume terms are 
taken into account while computing the full ChPT corrections to fx relevant to the lattice 
computation of this quantity with staggered quarks. 

In fig. 2, we illustrate the finite volume effects as predicted by the above formulae. Since 
the mass of the strange quark is directly simulated on the lattice, we keep it fixed. As 
for the light quark we define, r = m^/mP^^*^, which we then vary as r G [r„/rf,l], where 
rujd = {rriu + md)/2ms = 0.04 [20]. We also use the Gell-Mann-Oakes-Renner (GMOR) 
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and Gell-Mann-Okubo formulae, namely, 



ml = 2Bomsr , m'j^ = 2Boms^— , = 25om,^— . (36) 

The illustration in fig. 2 is made for realistic volumes currently used in lattice simulations: 
L G [1 fm, 2.5 fm]. In obtaining the quenched curves, we assume mo = 0.6 GeV and 
ao = 0.05. The curves are insensitive to the value of a^. On the other hand, the effect 
shown in fig. 2 becomes pronounced if mo is increased to mo = 0.65 GeV or mo = 0.7 GeV, 
values sometimes also quoted in the literature. 

From fig. 2 we see that the finite volume effects become more pronounced as the light 
quark gets closer to the physical u/d quark mass. In particular, they result in shifting the 
quenched fx toward a larger value, wheras the shift of fx in full (unquenched) QCD is 
opposite, i.e., the finite volume effects lower the value of fx- The partially unquenched 
cases lie between the two. We see that the quenched chiral logs become dominant as soon 
as the mass of the valence quark becomes lighter than the sea-quark mass. We did not plot 
the case when the valence and sea quarks are degenerate since such a curve is very close 
to the full ChPT case. It is worth noticing that in the region of r < 0.2 the finite volume 
effects for the partially quenched fx, with rsca = 1, are larger than for the quenched case. 



5.3 Finite volume corrections to Bk 

We proceed in a completely analogous way as in the last subsection and define 
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The corresponding ChPT expressions read 
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Fi gure 3: The finite volume corrections to Bk in full, partially quenched and quenched theory, 
eqs. (38), (39), (40), respectively. 
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The illustration, similar to the one discussed in the previous subsection, is provided in 
fig. 3. We observe that, as in the case of fx, the finite volume effects become pronounced 
as the light valence quark is lowered toward the physical u/d quark mass. Moreover, they 
have a tendency to enhance the nonlinearities which would otherwise be attributed to the 
physical chiral logarithms. 



^By "physical chiral logarithms" , we mean the chiral logarithmic behavior in infinite volume. 
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5.4 Asymptotic L ^ cx) formulae for fx and Bk 

As discussed at the beginning of this section, the difference of the one-loop ChPT formulae 
for the pion mass in the finite and infinite volume (i) reproduces Liischer's general formula 
at LO, and (ii) allows one to estimate the size of the corrections suppressed in Liischer's 
formula. Similar asymptotic formulae for fx and Bx can be obtained from our full ChPT 
expressions for Afx/ fx and ABk/Bk, i.e., eqs. (33) and (38), respectively. By using the 
asymptotic form (29), we get 

4VT; (27rm.L)3/2 ' ^ 2 m\ [T J (27rm.L)3/2 " ^ > 

To check how good an approximation these formulae are to the complete ones, given in 
eqs. (33) and (38), we made a numerical comparison of the two, and conclude that for 
volumes larger than (2 fm)^ and masses r > 1/4, eq. (41) is an excellent approximation. 
Otherwise, i.e., in the region in which the finite volume effects become important, the 
asymptotic forms (41) become inadequate and eqs. (33) and (38) should be used. We note, 
in passing, that formulae similar to those in eq. (41), but in the quenched case, were reported 
in refs. [5, 6]. 



6 Summary 

In this work we computed the one-loop chiral corrections to the decay constant fx and 
to the bag parameter Bx in all three versions of ChPT, i.e., full, quenched, and partially 
quenched. After working out the formulae in both infinite and finite volumes, we were able 
to discuss the impact of the finite volume effects on the chiral behavior of fx and Bx- We 
show that in most situations the physical chiral logarithms are completely drowned in the 
finite volume artefacts. In other words, unambiguously disentangling the physical chiral 
logarithms from the finite volume lattice artefacts does not appear to be feasible unless very 
large volumes are used. 

We also discussed the shift of fx and Bx induced by the finiteness of the volume. In our 
discussion we fix the strange quark in the kaon to its physical mass m^, whereas the mass 
of the light quark is varied between ms/25 and mg. This mimics the current lattice QCD 
studies in which the strange quark is directly accessed on the lattice while the accessible 
light quarks have r = mg/m^ G (0.5, 1), so that an extrapolation to the physical Vu/d = 0.04 
is necessary. 

The results of our calculation indicate that for r > 0.25 the finite volume effects are very 
small as long as L > 2 fm. In that region we provide a simple asymptotic formula which is 
an accurate approximation of the full ChPT expressions. 

From our formula it is also obvious that the finite volume corrections to fx and Bx are 
different in quenched and partially quenched QCD from those obtained in full QCD. There- 
fore, if in practical numerical simulations one wants to see only the effects of (un) quenching, 
the finite volume effects must be kept under control. 



14 



Acknowledgment 



It is a pleasure to thank Gilberto Colangelo, Stephan Diirr, Vittorio Lubicz and Guido 
Martinelli for discussions and comments on the manuscript. Partial support by the EC's 
contract HPRN-CT-2000-00145 "Hadron Phenomenology from Lattice QCD" is kindly ac- 
knowledged. 



References 

[1] C. Aubin et al. [MILC Collaboration], hep-lat/0309088. 
[2] M. Battagha et a/., hep-ph/0304132. 

[3] J. Gasser and H. Leutwyler, Nucl. Phys. B 250 (1985) 465. 

[4] J. F. Donoghue, E. Golowich and B. R. Holstein, ''Dynamics Of The Standard Model", 
Cambridge Monogr. Part. Phys. Nucl. Phys. Cosmol. 2 (1992) 1. H. Leutwyler, 
"Chiral dynamics", hep-ph/0008124; A. Pich, Rep. Prog. Phys. 58 (1995) 563 
[hep-ph/9502366] ; G. Ecker, "Strong interactions of light flavours" , hep-ph/0011026; 
U. G. Meissner, Rep. Prog. Phys. 56 (1993) 903 [hep-ph/9302247] ; G. Colangelo, 
G. Isidori, "An introduction to ChPT", hep-ph/0101264. 

[5] C. W. Bernard and M. F. Golterman, Phys. Rev. D 46 (1992) 853 [hep-lat/9204007] . 

[6] S. R. Sharpe, Phys. Rev. D 46 (1992) 3146 [hep-lat/9205020] . 

[7] C. W. Bernard and M. F. Golterman, Phys. Rev. D 49 (1994) 486 [hep-lat/9306005] . 

[8] S. R. Sharpe and N. Shoresh, Phys. Rev. D 64 (2001) 114510 [hep-lat/0108003] . 

[9] S. R. Sharpe, Phys. Rev. D 56 (1997) 7052 [Erratum-ibid. D 62 (2000) 099901] 
[hep-lat/9707018]. 

[10] M. F. Goherman and K. C. Leung, Phys. Rev. D 57 (1998) 5703 [hep-lat/9711033] . 

[11] J. F. Donoghue, E. Golowich and B. R. Holstein, Phys. Lett. B 119 (1982) 412. 

[12] J. Bijnens and J. Prades, JHEP 0001 (2000) 002 [hep-ph/9911392] ; S. Peris and 
E. de Rafael, Phys. Lett. B 490 (2000) 213 [hep-ph/0006146] . 

[13] M. F. Goherman and K. C. Leung, Phys. Rev. D 56 (1997) 2950 [hep-lat/9702015] . 

[14] I.S. Gradshteyn, I.M. Ryzhik, Table of Integrals, Series, and Products, 5th ed. San 
Diego, CA: Academic, 1994 (see page 921). 



15 



[15] M.J. Lighthill, Introduction to Fourier analysis and generalised functions, Cambridge 
Univ. Press, 1958 (see page 69). 

[16] G. Colangelo and S. Durr, hep-lat/0311023. 

[17] M. Luscher, Commun. Math. Phys. 104 (1986) 177. 

[18] J. Gasser and H. Leutwyler, Phys. Lett. B 184 (1987) 83. 

[19] G. Aubin and G. Bernard, Phys. Rev. D 68 (2003) 074011 [hep-lat/0306026] . 

[20] H. Leutwyler, Phys. Lett. B 378 (1996) 313 [hep-ph/9602366] . 



16 



